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Abstract. We consider a database separated into blocks. Blocks containing target 
items are called target blocks. Blocks without target items are called non-target blocks. 
We consider a case, when each target block has the same number of target items. 
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algorithm is based on Grover-Radhakrishnan algorithm of partial search. We minimize 
the number of queries to the oracle. We analyze the algorithm for blocks of large size 
and also shows how to use the algorithm for finite blocks. 
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1. Introduction 

Database search has many applications. Search algorithm enters as a subroutine in 
many important algorithms. Grover discovered a quantum algorithm, which searchs 
faster than classical [lj. It starts searching from a uniform superposition of all input 
states in the database Grover algorithm was proven to be optimal E]. If a 
database has one target item Grover algorithm can find it in 

flail = N - 00 C 1 ) 

queries [iterations] . If the database has z target items and we want to find one of them, 
Grover algorithm can find it in 

7T [N 

flull=4V7' N ^°° ( 2 ) 

queries to the oracle, see [3]. This is a full search. 

Partial search considers the following problem: a database is separated into K 
blocks, of a size 

A user wants to find one target block, which contains the target item. Grover and 
Radhakrishnan discovered a quantum algorithm for partial search 0j. A general 
structure of the algorithm is 

• Step 1. Global queries: j\ standard Grover iterations, see (j3J). 

• Step 2. Simultaneous local searches in each block: j 2 local Grover iterations, see 
©• 

• Step 3. One reflection over the average amplitude in the whole database, see (jOJ. 

Partial search also starts from a uniform superposition of all basis states (0). For many 
large blocks an improved version of partial search algorithm was suggested in [Hj. The 
number of queries to the oracle for several blocks of large size was minimized in jU]. This 
is GRK algorithm which we shall use in this paper. Other partial search algorithms were 
considered in [Zj. GRK algorithm was formulated in terms of group theory in [13J, this 
provides a powerful mathematical tool for proof of optimality. Partial search was first 
discovered while study a list matching, see 8J. The GRK algorithm was formulated for 
blocks of finite size in Let us emphasize that only one target block with one target 
item was considered in |IJ 03 EH E] • To see the universal features of the algorithm we 
consider large blocks b — > oo. As the size of the blocks increases the numbers of queries 
during each Step ji and j 2 are controlled by parameters r\ and a, see fl4"Uj) . We choose 
these parameters to minimize the total number of queries ji + j 2 , see (|47j) . 

In this paper we generalize GRK algorithm to the following. The database has N 
items, it is separated into K blocks. Each block has b items, so JV = Kb. Some of the 
blocks are target blocks, and we denote the number of target blocks by t. Each target 
block has several target items, and we denote the number of target items in a target 
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block by r (same for each target block). Total number of target items in the whole 
database is z = tr. We shall denote a target item by \m). A collection of all target 
items in the database is A = {|m)}, a collection of all target items in a target block is 
a = {\m)}. Non-target items are \x), a collection of all non-target items in the database 
is X = {\x)}, a collection of all non-target items in a block is x = So the whole 

database is D = A U X. Non-target block B = x, but target block B = x U a. We 
optimize Grover-Radhakrishnan algorithm to find one of the target blocks. 

This paper consists of two parts. In the first part we consider the limit of large 
blocks b — > oo. We optimize Grover Radhakrishnan algorithm to find one of the target 
blocks. We compare our result with GRK algorithm jHj and show that the optimal 
version of partial search in a database with several target blocks can be obtained from 
GRK [one target block with one target item] by replacing b by b/r and K by K/t , see 
fl48H5()j) . In the second part we consider application of partial search to finite blocks. 
For this case, we find a sure success way. 



2. The Partial Search Algorithm 

2.1. Step 1: Global Queries 

The Grover algorithm is based on two operations. First one 

tr 

I t = I -2^2\m)(m\ (4) 

radA 

selectively inverts the amplitude of the states of target items |m), where J is the identical 
operator. Initial state of the database is a uniform superposition of all basis states as 

1 N ^ 

|si) = ~7^Y1 I s )' (si|si) = l- (5) 

* x=0 

Here N is the total number of items in the whole database. Second operation 

I sl = i-2\ Sl )( Sl \ (6) 
selectively inverts the amplitude of the uniform superposition. Consider a vector 

N-l 

\v) = ^2a x \x). (7) 

x=0 

The operator — I S1 inverts the coefficients about the average as 

JV-l . N-l 

- I si \v) = y^ag|x), & x = 2a-a x , a = —^2a x . (8) 

x=0 x=0 

Finally, Grover iteration is a unitary operator as 

Gi = -I Sl I t . (9) 
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From now, we shall call it a global iteration. We have to apply it several times G 3 ^ . To 
describe global query more precisely we shall need an angle Q\ defined by 

sin 2 ^ = |. (10) 

We assume that initially our database is in uniform superposition (j3J) as 

During the first step we use only standard Grover iterations ©• We shall use 
eigenvectors of G\ as 

G x \il>t) = A±|^>, \f = exp[±2i0!], (11) 

They were found in The first Step of the partial search consists of j± standard 
Grover iterations applied to the whole database. After j\ queries the wave-function of 
the database is 

"M-.; Sin((2j ; + 1)h) ±\r-) , "f\r) (12) 

It also can be found in EI. 



2.2. Step 2: Local Queries 

During the second Step local searches are made in each block separately in parallel. We 
can write it as a direct sum with respect to all blocks 

G 2 = ®G 2 . (13) 

In each block local search is a Grover iteration as 

G 2 = -I s J t . (14) 



Here 



and I S2 is 



2^2\m)(m\ (15) 



I S2 = I -2\s 2 )(s 2 \, \s 2 ) = -t= <s 2 |s2> = 1- (16) 

one block 

Physical implementation of local search is similar to implementation the standard Grover 
search. Local search is the Grover search in a block. Remember that the number of 
items in each block is b = N/ K. First let us consider the wave-function of a non-target 
block [a projection of the wave function of the database to this block] as 

\nB) = a nt 2j \ x ) = a m Vb \s 2 ). (17) 

one block 

We can read the coefficient from (|12jl as 
cos ((2j! + 1)6,) 

- = VN -tr ■ (18) 
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The wave-function of a non-target block does not change during the second step as 

I t \nB) = \nB) = -I S2 \nB) = G 2 \nB). (19) 

Now let us consider the target block. After Step 1 its wave function can be obtained 
from (fT2|) as 



\B) 



sin ((2^ + 1)00, , , cos((2j' 1 + l)0 1 )Vfc 3 r 

—w + — 



Vt ""' ' y/N-tT 

Here |u) is the normalized sum of target items as 

I T 



ntt). 



and | ntt) is a normalized sum of all non-target items in a target block as 

b-T 



\ntt) 



1 



T 



(ntt\ntt) = 1. 



target block 



Note that we can express the state (jlfij) as 



\s 2 ) = ^-r \fi) + 
We shall need an angle 2 as 
sin 2 2 = l. 



6-r 



| ntt). 



(20) 
(21) 

(22) 
(23) 



(24) 



Eigenvectors of G 2 are constructed similar to eigenvectors of G\, see (fTTj) . as 

G 2 |^> = XtWt), A 2 = exp[±2i0 2 ], 

lV>a> =^=\ti)±^=\ntt). 
We can resolve these equations as 

\ntt) = -L (_|^+) + IV;-)) , |„) = -L (|^+) + |^-» . 

Now we can express the wave function of the target block in terms of the eigenfunctions 
of G 2 as 

^((2^ + 1)0!) 



\B) 



2t 

sin ((2^ + 1)0!; 



— i 



+ i 



b-T COS ((2^ + 1)0!) 



N-tr 



V2 



b-T cos ((2ji + l)^; 



'2t \ N — tr 

After j 2 iterations the target block will be 
G J 2 2 \B) = 

■"sin ((2^ + 1)00 



v/2 



l^->. 



(25) 



>/2 

,— 2102^2 



sin ((2^ + 1)0! 



— z 



COS ((2^ + 1)0! 



' W^—^cos ((2^ + 1)0!) 



m + 



(26) 
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Now we can use (J24j) to express the wave functions of the target block in the following 
form 

\B 2 ) = Gf\B) = at\fj) + a ntt \ntt). (27) 

and to calculate the amplitudes 

cos(2j 2 2 ) sin ((2.7! + 1)00 
a t = -= h 

Vi 

h ~ T sm(2j 2 9 2 ) cos ((271 + 1)00, (28) 



a>ntt — — 



N-tT 

sin(2j 2 2 ) sin ((2^ + 1)0! 



hl cos(2j 2 2 ) cos ((2jx + 1)00- (29) 



N-tr 
Later we shall also use 

h\B 2 ) = -atlfj) + a ntt \ntt) (30) 
2. 3. Step 3: Location of a Target Block 

The last Step consists of the following operation. We start by the application of the 
operator 

-/ ai = -J + 2| Sl )( Sl |. (31) 

[see (|SJ)] to the wave function of the whole database. This is a global inversion about 
the average a x — » 2a — a x , see (jEJ). 

The amplitude of an item which we want to vanish should be double of the average: 
a x = 2a. We want to annihilate the amplitudes of items in the non-target blocks 
a n t — ► 0. This means that the amplitudes introduced in the previous subsection should 
satisfy a nt = 2a as 

2 



ant ~ N L 



b(K - t)a nt + ty/ra t + Wb - ra ntt (32) 
see (J2TJ). We can rewrite this as 



K 



b y-— + tj a nt = t\Jra t + Wb - ra ntt . (33) 
Substituting here the expressions (fTHj). P9*jl we obtain the cancellation equation as 

trcos(2j 2 2 ) sin ((2^ + 1)00 + 
^/r^^sin(2j 2 2 )cos((2j 1 + 1)00 - 



y/t(b - r) sin(2j 2 2 ) sin {(2 3l + 1)00 + 

t h r^— cos(2j 2 2 ) COS ((2j! + 1)00 . (34) 
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This equation guarantees that the amplitude of each item in each non-target block 
vanishes as 

-I Sl \nB) = 0, (35) 
- I Sl \B 2 ) = -(a* - Vra nt )\fi) - (a ntt - \/b - ra nt )\ntt). 

Here we used (J27|) . (jHJ) and a nt = 2a. 

Now we can do a measurement. In the simplest case N = 2 n and K = 2 k , we label 
blocks by qubitsfitems inside of a block are labeled by n — k qubits]. We measure only 
k block qubits and find a target block. 

Let us make an extra query to the oracle to put a target block into a canonical 
form as 

\B 3 ) = —I t I sl \B 2 ) = sino;|/i) + cosu\ntt) (36) 
= (a t - \fra nt )\ij) - (a nU - \Jb- ra nt )\ntt). (37) 

In later, we shall calculate the angle u, see (}53|) . 

Meanwhile, we can also apply last two operations in inverse order as follows. 



\B 3 ) = -I sl I t \B 2 ) = (a t + \fra nt )\[i) - (a ntt - Vb - Ta nt )\ntt) (38) 

Now to vanish amplitudes of items in non-target blocks (}35|) we will have to change the 
sign of a t in equations (|3~2j) and (JBlIj) . compare (|^7|) and (JBUj) . Also the sign of first two 
terms in the right hand side of eq (J3~2j) . 

2.4- The Limit of Large Blocks 

To see universal features we consider the limit when each block is very large b — > oo, 
this makes the total number of items in the whole database also large iV = Kb — > oo. 
The expressions for angles (fTUjl . (j2SI) simplifies as 

The numbers of global and local iterations scales as, respectively 

ji = \ h = aA (40) 

similar to t = 1 , r = 1 case |U El E| • Here r] and a are parameters of order of 1 [they 
have a limit as b — > oo]. We shall minimize the total number of queries and find the 
best values of these parameters as a functions of number of blocks K, see This will 
give the optimal distribution of number of queries between Steps 1 and 2 of the partial 
search algorithm. We shall need the arguments of trigonometric functions in (|34jl as 

(2^ + 1)0! = | -2W^, 2j 2 e 2 = 2ay/¥. (41) 
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Let us leave only leading terms [of order a/6] in the equation as b — > oo 



v 7 ^ t \ . f2rjVtr 



sin 



if/ V Vif 



/— , f27]Vtr\ FT . „ / 

= — v ot sin(2aV' ir ) cos — + tW — cos(2av' r ) sin 

\ VK J V K v VV^ 

Note that at does not contribute in this limit, so this equation is the same for both cases 
of preparation of the final state of a target block and (JHHj) . 
We can simplify the equation to 

'VK t + tcos(2*VT)\ si j2WT\ 



2 y/K y/K J V VK 

n ■ , n /-> (2r)y/tr\ 
V t sm(2a\/Tj cos I ^_ 1 . 



We can express 77 as a function of a as 



(2^\ = ^si^Vr) 
\ VK J K-4tsin 2 (av^) 
see (HOJ). Notice that after we introduce rescaled variables 
K 

K = —, a = ay/r, r] = rjy/r, (44) 
the equation (J4*H|) will coincide with corresponding equation from [S| 



, 2fi \ 2VKwo.(2ol) 
A — 4 sin (a 



£.5. Optimization 

The most efficient algorithm makes the least number of queries to the oracle as 

jl+h + l^^^+(a-r ] )Vb^^^+^(a-fj). (46) 
To optimize the algorithm we have to minimize 

a-fj, 

with constraint (|45jl . The optimal values of a and 77 can be copied from [E]. The 
minimum number of queries can be achieved when 

2f}(K) .-, K-2 
tan /v ^ ; = , cos 25(10 = — = . 47 

This describes the optimal distribution of queries between Step 1 and Step 2, see (pfUj) . 
We shall consider these formulae as an expression of fj and a as functions of K. 
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Let us summarize the result. The optimal number of global queries at Step 1 is 

*-{W!-«(t)}£ 

Here the function 



depends only on the ratio Kj t. The number of local queries at Step 2 is 



. K\ 1 / K-2t , 
a — =7: arccos I — — . (51) 



where 

J ,_ „ 

After Step 3 amplitudes of all items in all non-target blocks vanish. 

Let us discuss the final wave function of the target block Note that a nt does 
not contribute in the coefficient in front of \fj) as b — ► oo see fjTKj) . so wave functions 
(j3T)|) and (J3~H|) coincide as 

|£>3) = \B 3 ) = sinc<j|/i) + cosc<j|ntt). (52) 

At the minimum (J47)) we get 

u = a (j J . (53) 



3. Sure success partial search 

5.1. Motivation and Idea 

Until now, we have focused on how to find a target block in limit of large blocks. The 
algorithm should be modified for case of finite blocks.. Hence we should find a way for 
a realistic case when b is finite. In this section, we explain a sure success way for this 
case. 

Briefly, we will not change the total number of queries to the oracle ji +J2 + 1, we 
will only change Step 3. In Step 3, we shall modify both operators 7 S1 and I t by phases 
in such a way that the algorithm will find a target block with 100% probability. 



3.2. Sure success way 

A sure success way for the case of single target item/single target block partial search 
was already proposed in [12] ■ Therefore we can extend this sure success way for our 
problem, partial search with multiple target items/multiple target blocks. As in [7] IT2"]. 
we note that the entire action of the partial search can be compactly described by a 
3-dimensional subspace spanned by the vectors: \M) the normalized state of all target 
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items; \NTT) the normalized state of all non-target items in all target blocks; and \u) 
the normalized state for all other items. Three bases are respectively 

E |^) = -^^|m), 
target blocks meA 

t „ t b-r 



\M) 
\NTT) 



1 

7t 



i 



E 



\ntt) 



\u) 



target blocks 
1 



E Ei-)' 



target blocks xex 



b(K-t) 

E 



;r) 



items in non-target blocks 
Using these orthonormal basis, the initial state |si) of (J5J) may be written 

\si) = sin 7 sin 9 2 \ M) + sin 7 cos 2 1 NTT) + cos 7 1 u) , 



(54) 
(55) 
(56) 

(57) 



where sin 2 7 = t/K , sin 2 6»i = tr/iV, and sin 82 = r/b . Note that these three bases 
states are generalized from |12j . From this, all subsequent flows are the same with the 
result in |T2j . 

Based on the three bases, (7f operator is represented as 



n 



TM h T, 



where 



/ 1 







q cos 82 sin 7 
cos 81 

^ 

COS 8± 



cos 7 
cos 81 
cos (?2 sin 7 
cos 8\ 



(5c 



(59) 



/ 



and 



M, 



.1 \ 



I cos{2j x e x ) sin^M) 
-sin(2ji6'i) cos(2j'i^i' 
\ 



\ 




(60) 



Lr 2 



Also, based on the same three bases, G J 2 operator is represented as 

/ cos(2j 2 2 ) sin(2j 2 2 ) \ 
-sin(2j 2 2 ) cos(2j 2 2 ) 
1 

The intermediate state after j\ global Grover iterations is given by 

/ cos^i (s g m + c g cos 6*i sin#i) \ 



V 



(61) 



/ 



cos 2 6*i 



(62) 



/ 



cos 6 2 sin 7 (c 9 m — s g cos 6\ sin #1) 
1 cos 7 (c 9 m — s 9 cos 6\ sin #1 ) 

where c 9 = cos(2ji#i), s 9 = sin(2ji6 l i) and m = cos 2 2 sin 7 + cos 2 7. For a 
more compact presentation of Gf\si), let k = s 9 m + c 5 cos #1 sin Q\ and / = c 9 m — 
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<%\*i) 



cos 2 Q\ 



\ 



cos 6*2 sin 7 (I) 

COSj(l) j 



(63) 



The next intermediate state after j% global and j% local Grover iterations is given by 
GfGf| Sl ) 



cos 2 9\ 



/ 


C\ cos 6i (k) + si cos 9 2 sin 7 (Z) 






/ 


a 


\ 




—si cos #1 (/c) + c\ cos # 2 sin 7 (I) 








b 




V 


cos 7 (7) 






V 


c 


J 



(64) 



where q = cos(2j 2 #2) and si = sin(2j 2 6 l 2). 

The final global Grover operator iteration, —I s I t , is modified with two phases as in 
the exact Grover search [Oj 

^final = _ [/ 



(I - e 2ie )\ Sl )( Sl \] x [I — (I — e**~*>)\M){M\] 



(65) 



Translating this into the three basis states supporting the entire computation we 
obtain 



^final 

V 



v 



A4>-6) 



1 [1 — p sin 2 7 sin 2 # 2 ] P sin 2 7 sin # 2 cos 9 2 
p sin 2 7 sin 62 cos #2 p sin 2 7 cos 2 #2 — 1 



3 i(0-0) 



Pf 



pg 



pf 
pg 

p cos 2 7 — 1 j 



(66) 



where p = 1 — e , / = sin 7 sin # 2 cos 7, and (7 = sin 7 cos 7 cos 62- 

Finally then, our aim of a sure success partial search will be achieved if we can 
find two phases, 9 and <ft, for the above final global Grover operator which satisfies the 
condition 

Kx|G^G?Gf|si)| = 0. (67) 

The relevant phase condition then reduces to 

ajtt-Vpf + bpg + c[p cos 2 7 - 1] = , (68) 

where a, b and c are defined in (fMj) . 

The phase condition may then be rewritten as 

e i{<f>-e) px + py + 2z = ; (69) 

where 



x = af, y = bg + ccos 2 7, z = -- . 

The real and imaginary parts of ()69)1 may be simplified to give 

V - - z 



(70) 



sin0= sin# : — -, 

x x sin t) 

cos 4> = — — cos 6 . 

x 



(71) 
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Finally, combining these two equations together, we may eliminate to yield 



z 



2 



sin 2 6 = — , (72) 

or — y z — lyz 

which to have a solution must satisfy 

x 2 > (y + zf . (73) 

There will then be a solution for <f> provided the right-hand-sides of (J7TJ) are bounded 
in absolute value by unity. 



4. Summary 

In the paper we described a quantum algorithm for partial search in a database with 
several target items. A database of N items separated into K blocks of b items each. 
The partial search algorithm finds a target block. We found a simple dependence on 
number of target blocks and on number of target items in each block. In the first half 
of the paper, we studied the algorithm for large blocks to understand universal features. 
In the second half of the paper we explained how to apply the algorithm in a case of 
finite-sized blocks for realistic applications. 
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